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Non-universal dynamics of dimer growing interfaces 
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A finite temperature version of body-centered sofid-on-sofid growth models involving attachment 
and detachment of dimers is discussed in 1+1 dimensions. The dynamic exponent of the growing 
interface is studied numerically via the spectrum gap of the underlying evolution operator. The 
finite size scaling of the latter is found to be affected by a standard surface tension term on which 
the growth rates depend. This non-universal aspect is also corroborated by the growth behavior 
observed in large scale simulations. By contrast, the roughening exponent remains robust over wide 
temperature ranges. 
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^ • I. INTRODUCTION 

(N ■ 

In studying statistical aspects of non-equilibrium surfaces the onset of scaling regimes at both large time and length 
(~| . scales has enabled one to characterize a vast body of growth processes in terms of universality classes jl| . In analogy 
O " to equilibrium phase transitions, there is consensus in that the late evolution stages of these processes are controlled 
^ by a set of scaling exponents stemming ultimately from the symmetries and conservations laws of the underlying 
growth rules. A basic quantity of interest investigated extensively in this context concerns the roughness or surface 

,^ width W{L,t) developed by growth fluctuations at a given time t when starting from an initially flat substrate of 

• typical length L. Based on a wide range of theoretical and numerical studies it can be argued that W scales as [l|, ^ 



W{L,t) = L<f{t/L^), (1) 



with a universal scaling function behaving as f{x) ~ x^^^ for a; -C 1 , whereas for x ^ 1 it remains constant. 
Consequently, for t ^ L'' the width saturates as while growing as t^^^ in the thermodynamic limit. The roughening 
'"^ ■ exponent C measures the stationary dependence of the surface width on the typical substrate size while z, frequently 
^ I referred to as the dynamic exponent, gives the fundamental scaling between length and time. 

O , In this work we focus on rather unusual scaling properties studied in recent years both in one 0, 0, Q smd two 
I ^ I ■ dimensional interfaces Q by means of discrete models of surface growth. For simplicity, here we consider a body- 
centered solid-on-solid (BCSOS) version of these [11 , limiting height differences between neighbors to ±1 rather than 
CN ■ to 0, ±1, as in restricted SOS reahzations [1, Our basic kinetic steps (depicted schematically in Fig.[T]), involve 
^ ' adsorption and desorption, possibly after recombination, of dimers only. Attempts of desorption can take place 
whether or not the selected pair of adjacent monomers arrived together. Therefore, the rule for evaporation allows 
for reconstitution of dimers, a crucial feature, so their identity is not maintained during the growth process, 
^y!^ It is important to note that throughout the stochastic evolution the parity of the number of monomers (eventually 

isolated) is conserved at each height level of the surface. The implications of this 'evenness' non local constraint 
on the scaling exponents are far reaching Q, and in the one dimensional {Id) case have been analyzed in terms 
of even visiting random walks. In this latter representation, interface configurations are thought of as Brownian 
^ paths whose locations (i.e. height levels) are visited an even number of times before the walk terminates on a 
!• • given time interval (here playing the role of the substrate length). In marked contrast to normal random walks, the 
. 5^ I evenness constraint introduces highly correlated movements giving rise to an anomalous (sub diffusive) mean square 
I displacement, which in the interface language means a saturated width scaling not as L^/^ but rather as L^^^ 0]. 
$H ' Also, the dynamic exponents z obtained numerically for these [1, 0, [E| and related globally constrained systems 
1 [§, are definitely different from those of usual monomer type interfaces, irrespective of the later evolving towards 
equilibrium or nonequilibrium stationary regimes (as exemplified respectively by the Edward- Wilkinson Q (EW) and 
Kardar-Parisi-Zhang ^ (KPZ) universality classes). These anomalous aspects of non local constraints should not be 
regarded as purely academic. In fact, in catalytic surface processes the interplay between the substrate geometry and 
the shape of the intervening objects does matter. In particular, dimers become relevant in the roughening dynamic 
of vicinal surfaces which only allows deposition and evaporation of diatomic molecules [q] . 

To our knowledge, there are no available phenomenological equations of growth (e.g. KPZ, EW), suitable to 
describe interface fluctuations arising from deposition-evaporation of composite particles. The manner in which non 
localities could be incorporated in that continuum limit is by far not clear. However, in probing the robustness of 
the dissociative dimer models referred to above some progress can be made by introducing a continuously tunable 
parameter without affecting neither the symmetries nor the conservation laws of their dynamics. Specifically, we 
consider a finite temperature extension of Id dimer growing interfaces, which for the ease of our numerical analysis 
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FIG. 1: The four relevant deposition-evaporation cases and their respective rates g , g~ for a one dimensional BCSOS dimer 
growing interface. Adding or removing a dimer (denoted by dotted lines) at columns hj,hj+2 can be viewed as flipping the 
spin-i quartet Sj, Sj+i, Sj+2, Sj+3. The involved rates depend on the change of surface tension ([2|) which is in turn determined 
by the spin (or slope) states of Sj^i and Sj+4 [see Eqs. (|3]) and Each situation is depicted schematically from (a) to (d). 



(Sec. Ill A), is here taken under detailed balance conditions. Following Ref. [11|], this is readily done through a 
standard surface tension term associated to configurational energy scales which discourage the development of strong 
fluctuations, besides those already prevented by both BCSOS and evenness constraints. In spite of these severe 
restrictions, they should not impede us to evaluate scaling exponents, as equilibrium surfaces at finite temperatures 
are always rough in Id. In the case of Id-monomer growing interfaces, the surface tension is not relevant to drive 
them out of their universality classes, though interest ing ly in 2d it can change the nonlinear term sign of the KPZ 
equation thus inducing rough-to-rough transitions pTi . Il2| | . Surprisingly, for dimers it will turn out that the combined 
effects of global constraints and surface tension entail significant changes in the z dynamic exponent, suggesting 
rather a non-universal temperature dependent value. In fact, non-universal aspects were already observed over wide 
temperature ranges in the monomer systems studied in Refs. [ill . [l2| , although they were ascribed to finite size effects 
which become particularly severe in nearing the equilibrium roughening transition temperature Tr p^ . However as 
mentioned above, in our Id case Tr is strictly zero and in practice all equilibrium correlation lengths can be fairly 
bounded as long as temperatures are not taken too low. In this sense, it is worth mentioning that steep variations 
of z will already appear within high temperature regimes (Sec. III). By contrast, the roughening exponent C remains 
robust as differences with respect to its vanishing tension limit (C ^ 1/3), will basically merge with our numerical 
errors. Since by definition C is a stationary index, it is just the nonequilibrium dynamics that is being strongly 
affected, as we shall see. 

Turning to methodological issues, as is known the dynamic scaling hypothesis referred to in Eq. ([T]) is usually put 
forward to determine both ( and z in a jointly way by simulating the growth dynamics over different substrate sizes. 
With the aim of obtaining an independent (separate) evaluation of these exponents, in addition to this standard 
procedure we will also exploit the known equivalence between BCSOS models and interacting gases of hard-core 
particles [l3] . Following the thread of ideas given in [l^ , we will recast the Metropolis operator that rules our growth 
simulations in terms of a quantum spin representation. This latter lends itself more readily for a finite size scaling 
analysis of the gap of the Metropolis operator which ultimately is related to the z dynamic exponent. On the one 
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hand this technique avoids the problem of deahng with long transient regimes though on the other is limited severely 
by the affordable substrate sizes. For now let us simply remark that already modest lengths are able to yield clear 
finite size trends over wide temperature ranges. To complement our approach, at low temperature regimes (where 
correlation lengths exceed the sizes reachable by exact diagonalization) , we will rely on numerical simulations of much 
larger systems which along with Eq. ([T]) will further support the non-universal picture. 

The layout of this work is organized as follows. In Sec. II we construct the quantum spin analogy of the standard 
Metropolis dynamic and briefly touch upon symmetries and conservation laws. By means of an ulterior non-unitary 
spin rotation, this results in a symmetric representation of the Metropolis operator. This simplifies considerably the 
subsequent numerical analysis of Sec. Ill in which the spectrum gap of this operator is obtained via standard recursive 
techniques p^ . The evaluation of dynamic exponents is then extended to low temperature regimes using standard 
Monte Carlo simulations. Finally, Sec. IV contains a summarizing discussion along with some remarks on extensions 
of this work. 



II. DYNAMIC AND REPRESENTATIONS 

As usual, the state of a solid-on-solid interface is represented by a set of single-valued functions {hj{t)} denoting 
height levels at positions j = 1, ... , L measured at a given time t from a reference substrate of length L. As mentioned 
earlier, to prevent arbitrary bulk fluctuations we impose BCSOS constraints on these heights which hereafter are 
taken to satisfy hj+i — hj — ±l,Vj, i, along with periodic boundary conditions (PBC). Growth or evaporation 
of the interface involves two particles (dimers) at the top of columns hj , which, to comply with the above 

restrictions, ought both to be local extrema of the evolving interface. More specifically, deposition (evaporation) 
events hj, hj^2 ~* hj + 2, + 2 {hj, ~^ hj — 2, ~ 2 ) can only occur at two consecutive local minima 
(maxima) of the heights set (see Fig.[T]). We stress that evaporation takes place regardless if these maxima were 
created together or not, so dimers can dissociate. Also note that the number of heights at a given level preserves its 
parity throughout. 

We want the transition rates of these processes to depend on the surface tension a referred to in Sec. I. In turn, any 
model that associates energies to height differences should provide a plausible description of cr. Due to our BCSOS 
choice {Ah — ±1), evidently the simplest form of a should assign an energy e > to each double facet exposed 
between columns hj and /ij+2- Since /ij+2 — hj = 0, ±2 , this may be studied by defining 

= f E I ^^^+2 -h,\, (2) 
j 

which simply counts the total number of double facets in a given interface configuration. Therefore, we can construct 
a standard Metropolis process at temperature T after introducing the following transition probability rates 

g^ihj , h,+2 h,±2, h,+2 ± 2) = min { e-^'"^^, l} , (3) 

where Aja is the change in surface tension upon depositing (5^) or removing {g^) a dimer at hj, hj+2 (henceforth, 
the Boltzmann constant is set equal to one). The four possible scenarios determining the values of these rates, 
which by construction obey detailed balance (that is g^{Aa/T)/g^{—Aa/T) = e~^'^^'^ ), are schematized in Fig.[T] 
Clearly, in the high temperature limit we recover the dimer model considered in whereas T — > 0"*" serves to favor 
smooth states over long transient regimes. 

It is often more practical to work in terms of slopes rather than with interface heights, so in what follows we will 
employ the known mapping between BCSOS and hard core particle dynamics [14.] . This correspondence is easily 
visualized in Fig. [1] which simply associates the height differences hj — hj^i = Sj to particles {sj = 1) or vacancies 
{sj = — 1). Conversely, the interface heights are obtained as hj — X]n<j ' modulo a constant level. In particular, 
in this picture the surface tension reduces to the Ising Hamiltonian 

3 

up to an irrelevant constant; whereas after some straightforward manipulations, the square deviation of the instan- 
taneous average height /i of a particular slope configuration |s), i.e. a given realization of the interface 'width' 
W^s) = T.j{hj -hf/L, has the form 

^^)-^^ + l2Y.'^L-^)'^'^- (5) 
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Creating (eliminating) a dimer now amounts to a backwards o«o» • o • o (forwards •o»o — > o»o») move 
of two particles at a time, whereas from Eq. ^ it can be readily checked that in this representation the corresponding 
rates ([3]) just involve the neighboring states of the flipped quartet Sj, Sj+i, Sj+2, Sj+3i namely 

gf = min | ^ " "^-^ \ 1 | , (6) 

with K = e/T being from now on our inverse temperature parameter. Note also that under PBC the interface 
'magnetization' Sj vanishes at all evolution stages. 

In addition, some constants of motion can be immediately identified. Evidently, the dimer dynamic exchanges two 
particles between four consecutive sites while changing the occupation of the involved next nearest neighbor locations 
by the same amount. If we think of these sites as being part of a four-partite lattice A = Ai-t-...-t-A4 {L/2 even), hence 
upon defining Sa = X^jeA '^J magnetization of sublattice Aq it is clear that the set of dependent quantities 

{ {—l)°'Sa — (—1)" Sa' , a, a' = 1, ...,4} (of which only three are independent), is left invariant throughout. From 
a more fundamental point of view these conservations arise ultimately from continuous symmetries borne by the 
Metropolis operator controlling the probabilities of our slope states, and towards which we now turn. 



A. The Metropolis operator 

As is known, the evolution operator of a Markovian process of the kind discussed so far can be constructed generically 
as [13 

{s'\M\s) ^ -R{s~^s') , sj^s', (7) 
{s\M\s) = J2 ^(^-^'), (8) 

s'^s 

where R{s — > s') denotes the transition rate at which configuration \s) evolves to |s') per unit time. At least formally, 
this enables one to derive all subsequent probability distributions \P(t)) = J^s^i^^^) \^) from the action of the 
evolution operator on a given initial state, that is \P(t)) = e~*^*|P(0)) In our problem, R{s — > s') = 1, e"^'' 

for all compatibles |s), |s'), and the specific form of M can be readily found by interpreting the slope configurations 
|s) = |si, ... , sl) as eigenstates of the z component, say, of Pauli matrices (Ti, ... , <7l assigned to each slope site. For 
instance, using spin-^ raising and lowering operators cr"*", , the operational analog of Eq. ([7]) will then read 



^ M,,. \s') (s| = - E ( at A + 9j A, ) , (9) 

s,s' . s^s' j 

where the adsorption (desorption) or double exchange operator Aj (Aj) acting on the j-quartet referred to above is 
simply 

^] = <^7+i<+2'^7+3- (10) 

Here, are thought of as diagonal operators in the representation and whose elements are identified with the 
rates of Eq. ([6|). In this regard, the ordering of application in ([9]) is immaterial. As for the diagonal elements of Eq. ([8]), 
needed for conservation of probability, they basically count the number of ways in which a given configuration \s) can 
evolve to different states \s') by fiipping an active quartet at a time. This can be properly tracked down in terms of 
number operators n — a~^(T~ = (1 + cr^)/2 and weighting each probed quartet with its corresponding rate (g^). The 
counterpart of Eq. ^ then becomes 

M,^, \s) (s| = g+ (1 - fij) fij+i (1 - nj+2) nj+3 (11) 

s j 

+ 9j (1 - "j+l) (1 - nj+3) , 



which along with Eq. ^ completes the form of our Metropolis operator. Taking into account the spin algebra 



[a^ , a, ] — Sij a-|, {af , cr }+ = 1, the former finally reduces to 



M 



j:{st4+SiA,){A]+A,-l). (12) 
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By construction M is a stochastic operator and therefore its ground state \^o) has vanishing eigenvalue and corre- 
sponds to the detailed balance solution of the problem, i.e. |\l/o) oc e~'5'^{=> |s). Instead, its left ground state {tp\ is 
an equally weighted linear combination of all reachable (s| (note that M is a non hermitian operator whose columns 
add up to zero). With the aid of this left state and starting from an initial probability distribution |P(0)), typical 
quantities of interest, such as the interface width W , are calculated as >Ve-*=^*|P(0)) Here, the 'width 

operator' VV is obtained by promoting the slopes of Eq. ^ to cr^ matrices. 

Despite the apparent simplicity of our Metropolis operator, the A's above can not be associated to elementary 
excitations of any kind and exact analytic treatments may seem unlikely. Nevertheless, Eq. p2p will permit some 
numerical progress on finite size systems after considering a simple transformation to be discussed later on in Sec. 
lie. Before that and for the sake of completeness, we pause to digress briefly about symmetries and conservation laws 
of M. 



B. Excursus: constants of motion 

Here we follow Refs. [3, [3 in closely related processes. Recalling that under a rotation by an angle 9 around the 
^-direction transform as e^'^cr^, we can therefore choose angles 9a for all spins in each sublattice Aq, such that 
Eq. (jlOp [and obviously pT|) ] is left invariant. Clearly, this is the case of 

j2{-ir9a = o. (13) 

a 

On the other hand, the infinitesimal generator of this transformation is = J2ai~^)°' , with ~ X^jeAo '^j- 

Since M = e''^/^Me~''^/^, then [M , S] = , and therefore S is preserved by M. But from the constraint ([TS]) it 
follows that S can be rewritten in terms of three independent angles, that is S" = J2a^a' — (—1)" S^,] 

from which one recovers the three conserved quantities identified before on more intuitive grounds. 

These continuous symmetries entail a number of invariant subspaces growing at most as L^, which however by 
no means exhaust all possibilities. For instance, it is straightforward to see that already the number of jammed 
configurations (i.e. states that can not evolve further), grows exponentially with the system size [18]. This unusual 
proliferation of invariant states should be the consequence of a much higher symmetry of M. Although its explicit 
operational form might be difficult to figure out, we can nevertheless follow Ref. [l9| and construct an exponential 
number of dynamically disjoint sectors, either jammed or unjammed, regardless of the value of k. To this end, one 
defines a reduction rule by looking at the occurrence of groups of active quartets in a given configuration \s). Each 
occurrence, if any, is deleted so the length of the remaining object is reduced in 4-bits per deletion. This procedure 
is applied recursively until one is left with a string that can not be further reduced, i.e. an irreducible string I{s}. In 
turn, the result is unique irrespective of the order of deletion. To mention only a few examples: I{\ TiT-iT-iT-iiT)} = iT 
(either in one or two steps); the flat interface or antiferro state yields a null string, whereas any jammed configuration 
is already an irreducible string of length L. The key issue to bear in mind is that two states \s) and \s') belong to the 
same M-subspace I{s} — IW} [19]. So, this non-local construct picks out both the length C ^ L — Ak and the 
sequence of the irreducible string's elements (their combinations growing exponentially in C), and ultimately defines 
the constant of motion under which the Af -dynamics take place. 

That being said, from now on we shall content ourselves with studying just the null string subspace selected by 
initially flat conditions, for the most part quite natural in the context of growing interfaces. Note also that the 
equilibrium properties of the Id-Hamiltonian ([¥]) are neither analytically simple to evaluate (e.g. Eq. ([51), not even 
for K = 0), as the ensemble of averaged states must be consistent not only with = (PBC) but with a vanishing 
irreducible string as well (totally unjammed conditions), which rules out an exponential number of states. 

C. Symmetric representation 

Returning to the discussion of Sec. HA, we may make some progress on the numerical analysis of M by performing 
a similarity transformation so as to map this operator into an hermitian matrix. This is feasible because detailed 
balance in rates ^ ensures the existence of a representation in which the evolution operator is self adjoint [l^. For 
this purpose, it suffices to consider a diagonal transformation alike the one discussed in Sec. II B but using pure 
imaginary angles instead. Specifically, we rotate each jth spin around the z-direction by a site dependent angle (field 
operator) 



= Y (^1-1 +'^1+1 ) 



(14) 
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by means of the nonunitary similarity transformation U = e "^^"^ . Under this rotation cj^ — > ^ (s~i + '^3+i)o-r^^ 
so it is simple to check that the double hopping operators of Eq. pHl) transform as 

U Aj U^^ = e"* - ''^+^)Aj . (15) 

This introduces new diagonal operators dj in the <t^ representation, that like the operators commute with Aj, Aj , 
and in terms of which the off diagonal part of M becomes symmetric. More specifically, this symmetrization is 
produced by defining 

dj = (""^1 " '"^'+^) g+ = e"* ('^^i - ''^+^) gj , (16) 

their diagonal elements being e~'5 ^^^-'^ ~ . As a result, Eq. Q is transformed into — dj{A'^- +Aj) while Eq. (fTTI) 
is left unchanged, so the rotated Metropolis operator H = UMU^-^ can be finally cast in the symmetric form 

^ = E (4 + ^^){97A] + at A, - 4) ■ (17) 

j 

Consequently, the time dependent probability distribution turns out to be a superposition of orthogonal eigenlevels 
■(Aa^o) with real eigenvalues A > of iJ, each having typical lifetimes 1/A. In particular, the ground state 1-00) has 
eigenvalue A = 0, and is just the transformed Boltzmann distribution j^'o) referred to above, i.e. \ipo) = C/|^o) « 
e^i'^*^) \s). Since left and right levels now coincide, it is thereby a simple matter to check that in the symmetric 
representation the dynamic of any diagonal observable, say the interface width Wl = X]i<j "i {L — j) erf (Jj (in 
turn invariant under U), can be written as 

W^{L,t) = W!g+ e-'-*(^ol>VLN'Aj(^o|J"(0)), (18) 

A^>0 

where W^^ is the saturation width reached at equilibrium, whereas |P'(0)) = U \P{0)) denotes the transformed initial 
distribution. From here we see that if the spectrum gap vanishes as 1/ , then a finite size scaling analysis of the first 
excited levels A^ will provide the dynamic z-exponent ruling over the late roughening stages referred to in Sec. I. 

It is worth pointing out that the discussion presented so far can be readily extended to include monomers (m = 1), 
trimers (m = 3), etc. , so long as the operators involved in Eq. (|17p are reinterpreted as 

m 
i=l 

gf = min{e±'^(^^+^"-"^-i),l} , (19) 

dj = g-f k,-l-S, + 2™| ^ 

In particular, for monomers with no surface tension (k = 0) the evolution operator reduces to the fully isotropic 
Heisenberg ferromagnet, thus recovering the usual EW dynamic exponent z = 2. Among other numerical aspects, in 
what follows we shall focus on the evaluation of this quantity at k > for both m = 1 and 2. 

III. NUMERICAL RESULTS 

The above ideas provide an alternative manner to evaluate dynamic exponents, independently of those obtained 
by the application of the dynamic scaling hypothesis ([T]). Thus, we first explore the consequences arising from the 
exact diagonalization of Eq. (jl7p in small systems, and then go on to corroborate them over larger length scales [via 
Eq.(IT])] using standard Monte Carlo simulations. In addition, these latter will complement the evaluation of z in low 
temperature regimes where, as we shall see, resorting to small lattice sizes might become inadequate. 

A. Scaling the gap 

To analyze our stochastic matrix we first obtained its null string basis using dimers on rings of sizes L = 4k. 
This was easily implemented by applying H (or alternatively, M) to either of the two antiferro states and keeping 
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proper track of the new generated configurations. By iterating this procedure with those new states for which H was 
not previously applied, the whole null-string subspace was finally expanded. The total number of states so found is 
of course independent of k, and as expected (see Sec. 11 B), grows slower than the monomer space dimensionality 
(l/2) °^ -^^^ M^ore specifically, this dimension seems to increase as « 1.6(7)^, at least for the sizes at hand, which in 
turn allowed us to explore rings of up to 32 sites [131 ■ 

Once having identified the null string configurations, we proceeded to evaluate exactly the low lying levels of H 
(in principle, just the first excited will do), via a recursion type Lanczos algorithm [16]. Starting that recursion 
from a random linear combination of null strings but chosen orthogonal to the Boltzmann-type distribution |^/;o) 
referred to above, we then obtained the finite size behavior of the dimer gap, i.e. of Ai. This is shown in Fig. [2] 
for several temperatures within the range < k < 1 indicating a gap decrease oc , however notice that the 
data do not fall into parallel straight lines. In an attempt to standardize this situation we used scalings of the form 
Ai(L) = Af^L^^^ + Bf^L~^^ , but a large amount of miccrtainty in both zi and Z2 raised doubts about the adequacy of 
such a procedure. Also, logarithmic corrections were attempted but no evidence supporting these latter were found. 
Thus, in principle we are led to suggest a plain power law decay although with a non-universal temperature dependent 
dynamic exponent z = z (k). In particular, z (0) ^ 2.6(1) is in fair agreement with the value obtained in Ref. [5] by 
standard simulations of dimer interfaces without surface tension. 
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FIG. 2: Finite size behavior of the first excited level Ai of the Metropolis operator (|17p . The inverse temperature k decreases 
from top to bottom, each symbol standing respectively for k — 1, |, ^, ^, 0. The dynamic exponents z are read off from the 
slopes of the fitting lines. 



A slightly improved estimation of z can be made by defining an effective dynamic exponent 

_ ln[Ai(£-4)/Ai(L)] 
^" ln[L/(L-4)] ' 

and then extrapolating to L — > oo for a given k. The results of this are exhibited in Fig. [31 which for comparison 
also displays the corresponding monomer dynamic exponents. These latter were derived using Eq. p9p for m — 1 
along with a similar numerical analysis but employing L = 2k {< 24) instead. Clearly, an EW behavior characterized 
by the Heisenberg exponent z = 2 mentioned a little earlier can be discerned in monomer interfaces, as opposed to 
dimer exponents which evidently are non-universal, at least if we are to judge by their ~ 20% variation between k « 
and 1. Although it is true that size effects increase monotonically our error margins with k, they are nevertheless 
fairly bounded within the range inspected (see also Fig. [2]). 

A measure of these size effects is provided by the equilibrium correlation length of the associated Ising antiferromag- 
net appearing in Eq. ([¥]). If this length becomes comparable to our available sizes, particularly at low temperatures, 
then the asymptotic dynamics will be distorted on approaching equilibrium as the average antiferromagnetic domain 
sizes, representing active regions in the interface, will be cut off by L. This situation is illustrated in Fig. [4] where 
we show the pair correlations C(r) = '^ji'^j ^j+r) evaluated in the ground state or equilibrium distribution of H 
[as stressed above, notice that analytic treatments are difficult even in this simpler case because of the null string 
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FIG. 3: Non-universal exponents of the dimer dynamics. Rhomboids indicate the results obtained from finite size scaling 
extrapolations of Metropolis gaps (see Fig. [2]). For comparison, the triangles exhibit these exponents in the monomer case 
(close to z « 2). Filled circles (joined by dotted lines) denote the z values arising from the dynamic scaling hypothesis ([T| 
applied to much larger systems at lower temperatures. 



constraint imposed on Eq. ([1]) ] . In between < k ^ 1 , traces of antiferromagnetic short range order are nearly 
smeared out and correlation lengths become small. However, above k ~ 2 they rapidly grow up and eventually get 
comparable to our maximum sizes, so precluding further analyses of the gap (which in fact comes out to be almost size 
independent). Thus, to complement the results obtained so far and check whether non- universal exponents actually 
extend down to low temperatures regimes, we finally turn to the dynamic scaling hypothesis and simulations using 
larger substrates. 
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FIG. 4: Pair correlations in the ground state of operator (|17|) for 32 spins. Filled circles {k — 0) are representative of high 
temperature regimes (which are quantitatively similar up to k ~ 1 ). Above k — 2 (open circles), the underlying correlation 
lengths become rapidly comparable to the lattice size. 



B. Simulations 



Following the Metropolis rules referred to in Sec. II, we evolved initially flat interfaces with L — 2^°, 2^^ and 
2^^ heights until reaching their stationary states. After a sequence of L update attempts at random locations, the 
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timescale was increased in one unit, i.e. t — > irrespective of these attempts being successful or not. Measurements 
of W{t) were carried out for k = 1.3, 2.5, 4 and were averaged typically over 10'* independent histories. In Fig. [5] we 
display one of the characteristic scaling curves obtained using Eq. ([T|) for k = 4 There, the data collapse was attained 
by setting roughening exponents C w 0.3(4) which are practically common to all temperatures studied (see also Fig. [7] 
below). By contrast, this is not the case of the z exponents which, in line with the results of Sec. Ill A, are severely 
altered by k. Although their precise values are blurred by our not too sensitive collapse conditions, nevertheless they 
do follow the non-universal trend already found with our gap analysis, as can be seen in Fig. [31 




FIG. 5: Dynamic scaling of the dimer interface width [Eq. ([1} ] taking k = 4. Sizes L = 2^^, 2^^ and 2^" are denoted respectively 
by triangles, circles and squares. The data collapse was attained upon setting ( ~ 0.3(4) and z 1.6(7). The dashed line is 
fitted with slope 2/3 = 2C/z. 

To corroborate further the validity of this claim, we also conducted simulations in much bigger scales measuring 
directly the growth exponent /3 = C/^- The reader's attention is now directed to Fig. [S] where the width evolution 
is contrasted at high and low temperature regimes in substrates of 10® sites. As expected, non-universal aspects 
show up: after averaging over ^ 40 histories, clearly two rather different /3 exponents emerge and hold for at least 
two decades. On the other hand, using the universal roughening exponent C ~ 1/3 already identified (see further 
estimations below), we thus obtain values of z consistent with those previously encountered in smaller systems. For 
comparison, the inset of Fig.[S] also shows the typical EW values of the corresponding monomer cases which, alike 
their dynamical exponents in Fig. [31 remain robust under wide temperature intervals. 

In respect of the roughening exponents, we finally considered the saturation or equilibrium widths Weq of a variety 
of substrate sizes subject to k = 0, 1, 2 and 4. The employed relaxation times cx L^'-'^\ range from 2 x 10® to lO'^ 
Monte Carlo steps for the largest cases and, as expected, decrease monotonically with k. Our results are displayed 
in Fig. [71 clearly suggesting a common value of C,. Due to the pair correlations involved in Weq [see Eq. (0], here 
size effects are also more noticeable at low temperatures. Nonetheless, a simple numerical fit of both amplitudes and 
slopes indicates that most of our data {L > 200) can be accounted for by the parametrization 

VFeg« 0.4(2) e-^'^L^, (21) 

with C « 0.3(2). This means that in equilibrium the interface becomes actually rough so long as InL ^ k is held in 
the thermodynamic limit. 

IV. CONCLUDING DISCUSSION 

To summarize, we have studied numerically the dynamics of dimer growing interfaces at finite temperatures using 
two independent procedures. The first one analyzes the spectrum gap of the evolution operator [Eq. ([TT]) ] by exact 
diagonalization of small systems, thus picking out dynamic exponents in a direct manner. Clear finite size trends were 
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FIG. 6; Growth of interface width for 10® sites at k = (upper curve), and k = 6. Dashed lines are fitted with slopes 2 (i. The 
inset displays the corresponding situations for monomer growing interfaces. 
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FIG. 7: Finite size behavior of the saturation width. The symbols stand in turn for k = (circles), k = \ (squares), k = 2 
(triangles) and, k = 4 (rhomboids). Solid lines display a common slope 2C,. As is shown in the inset, the amplitudes of these 
latter exhibit an exponential decay with slope —2^ [see Eq. (|2ip ]. 



obtained in the range < k < 1 (Fig. [21), and fairly bounded extrapolations were derived for z (Fig.[3|). Although 
the non-local symmetries (i.e. irreducible strings of Sec. II B) of our stochastic operators are unaffected by surface 
tensions, surprisingly the z exponents are non-universal, being dependent on k. Yet, a theoretical interpretation of 
such puzzling behavior remains quite open. This is in marked contrast with the dynamics of monomers interfaces, as 



their tensions do not take over neither the EW nor the KPZ universality classes, at least in 1 -I- 1 dimensions [11 [. 

Secondly, using the standard scaling hypothesis [Eq. ([T]), Fig. [5] we checked out these findings under lower temper- 
ature regimes where correlation lengths become larger than our maximum diagonalizable sizes (Fig. 3]). Despite the 
limited precision of this method for k > 0, our results confirmed the non-universal tendency observed in Sec. Ill A . In 
turn, measurements of growth exponents in much larger substrates (Fig.[6|) further validated the monotonic decrease 
of 

As for the roughening exponents (Fig. [7]), in all studied cases with k > the global constraint referred to in Sec. 
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I and further examined in terms of irreducible strings [18|, |19| , led to anomalous motions of a rather unconventional 
type (C ~ 1/3) J as compared to Levi flights and other restricted random paths i21(]. More specifically, they are 
consistent with those of even visiting random walks and not comprehensible in terms of EW or diffusive interfaces 
(e.g. monomers), which are definitely rougher. Under surface tension the range of correlations so introduced in the 
associated walk is finite, and therefore the scaling of its width must remain unchanged (consult Rcf. [21] ), though in line 
with Eq. (|21|) . its proportionality constant might depend on the precise form of these correlations. To endow further 
this robustness of C, it would be interesting to elucidate whether the analogy of non-interacting electrons moving in a 
random medium studied in Ref. [3] could be extended to the finite temperature interfaces (walks) investigated here. 

Other pending issues of interest concern starting the growth process from more general initial conditions (i.e. not in 
the null string sector), capable of modifying asymptotic regimes [sj, as well as considering biased dynamics (without 
detailed balance) such as those analyzed in Refs. [lll,[l3l- In principle, the first situation could also be studied with the 
methodology of Sec. II A; however for the second one the similarity transformation of Sec. II C is no longer useful and 
the unsymmetrization of the Lanczos recursion would be inevitable [IQ] . Finally, 2 + 1 dimensional generalizations of 
this study could shed light on the combined role that dimer dynamics and substrate geometry might have in catalytic 
processes. Whether or not non-universal aspects would also emerge there under surface tension, deserves further 
investigations. 
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